Search for: All records

Abstract This paper presents a fast and robust numerical method for reconstructing point-like sources in the time-harmonic Maxwell’s equations given Cauchy data at a fixed frequency. This is an electromagnetic inverse source problem with broad applications, such as antenna synthesis and design, medical imaging, and pollution source tracing. We introduce new imaging functions and a computational algorithm to determine the number of point sources, their locations, and associated moment vectors, even when these vectors have notably different magnitudes. The number of sources and locations are estimated using significant peaks of the imaging functions, and the moment vectors are computed via explicitly simple formulas. The theoretical analysis and stability of the imaging functions are investigated, where the main challenge lies in analyzing the behavior of the dot products between the columns of the imaginary part of the Green’s tensor and the unknown moment vectors. Additionally, we extend our method to reconstruct small-volume sources using an asymptotic expansion of their radiated electric field. We provide numerical examples in three dimensions to demonstrate the performance of our method. 

This paper presents a numerical method for solving the inverse problem of reconstructing the shape of periodic structures from scattering data. This inverse problem is motivated by applications in the nondestructive evaluation of photonic crystals. The numerical method belongs to the class of sampling methods that aim to construct an imaging function for the shape of the periodic structure using scattering data. By extending the results of Nguyen, Stahl, and Truong [Inverse Problems, 39:065013, 2023], we studied a simple imaging function that uses half the data required by the numerical method in the cited paper. Additionally, this imaging function is fast, simple to implement, and very robust against noise in the data. Both isotropic and anisotropic cases were investigated, and numerical examples were presented to demonstrate the performance of the numerical method. 

Abstract This paper is concerned with the inverse problem of determining the shape of penetrable periodic scatterers from scattered field data. We propose a sampling method with a novel indicator function for solving this inverse problem. This indicator function is very simple to implement and robust against noise in the data. The resolution and stability analysis of the indicator function is analyzed. Our numerical study shows that the proposed sampling method is more stable than the factorization method and more efficient than the direct or orthogonality sampling method in reconstructing periodic scatterers. 

This paper is concerned with the numerical solution to the direct and inverse electromagnetic scattering problem for bi-anisotropic periodic structures. The direct problem can be reformulated as an integro-di erential equation. We study the existence and uniqueness of solution to the latter equation and analyze a spectral Galerkin method to solve it. This spectral method is based on a periodization technique which allows us to avoid the evaluation of the quasiperiodic Green's tensor and to use the fast Fourier transform in the numerical implementation of the method. For the inverse problem, we study the orthogonality sampling method to reconstruct the periodic structures from scattering data generated by only two incident fields. The sampling method is fast, simple to implement, regularization free, and very robust against noise in the data. Numerical examples for both direct and inverse problems are presented to examine the efficiency of the numerical solvers. 

We consider the inverse problem of determining the geometry of penetrable objects from scattering data generated by one incident wave at a fixed frequency. We first study an orthogonality sampling type method which is fast, simple to implement, and robust against noise in the data. This sampling method has a new imaging functional that is applicable to data measured in near field or far field regions. The resolution analysis of the imaging functional is analyzed where the explicit decay rate of the functional is established. A connection with the orthogonality sampling method by Potthast is also studied. The sampling method is then combined with a deep neural network to solve the inverse scattering problem. This combined method can be understood as a network using the image computed by the sampling method for the first layer and followed by the U-net architecture for the rest of the layers. The fast computation and the knowledge from the results of the sampling method help speed up the training of the network. The combination leads to a significant improvement in the reconstruction results initially obtained by the sampling method. The combined method is also able to invert some limited aperture experimental data without any additional transfer training.